FUNCTION THEORY OF ANTILINEAR OPERATORS 
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Abstract. Unlike in complex linear operator theory, polynomials or, more generally, Laurent 
series in antilinear operators cannot be modelled with complex analysis. There exists a corresponding 
function space, though, surfacing in spectral mapping theorems. These spectral mapping theorems 
are inclusive in general. Equality can be established in the self-adjoint case. The arising functions 
are shown to possess a biradial character. It is shown that to any given set of Jacobi parameters 
corresponds a biradial measure yielding these parameters in an iterative orthogonalization process 
in this function space, once equipped with the corresponding structure. 
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1. Introduction. Classical complex analysis is ubiquitous in complex linear 
Hilbert space operator theory; see, e.g., [HI [TBI [T] and references therein. The inter- 
play between these two fields has found many forms of which the spectral mapping 
theorem is the most well-known and undoubtedly most widely used. In real linear 
operator theory, complex linear operators constitute one extreme while antilinear op- 
erators constitute the other one. There is strong evidence to expect that these two 
extremes turn out to be mathematically equally rich [T2j [131 [lOl [H] • In this paper, a 
function space structure is identified which has an analogous connection with antilin- 
ear operators as complex analysis has with complex linear operators. This function 
space surfaces in spectral mapping theorems for antilinear operators. An theory for 
it arises in the iterative construction of invariant subspaces for self-adjoint antilinear 
operators once the appearing Jacobi operators are connected with biradial measure 
spaces. 

For the spectral mapping theorems, suppose A is a bounded antilinear operator 
on a complex Hilbert space H. Antilinear means that A\ = \A for any complex 
number A. Then, for simplicity, take a polynomial p{\) = X]fc=o '^k/^^ ■ If x G if is an 
eigenvector of A, there holds 

p{A)x^p{\)x, (1.1) 

where p has transformed to p{\) = X]fe=o('^2fe + a2k+i^) \Xf'^ ■ More generally, by 
allowing more complicated analytic functions /, the transformed functions can be 
expressed as 

f{\)^u{\\\'') + v{\\r)\. 

with sufficiently regular complex valued functions u and v. These functions, denoted 
by C(r2), constitute a vector space over C carrying a natural notion of product. The 
spectral mapping for the spectrum a {A) then takes the form 

f{a{A))(Za{f{A)). (1.2) 
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The equality cannot be established in general. However, if A is additionally self- 
adjoint, the equality is shown to hold in (|1.2I) . In particular, when treated antilinearly, 
Hankel operators fit in this category in a natural way. 

For the theory for these functions, even if x is not an eigenvector of A, collecting 
all the possible vectors on the left-hand side of (II. 1|) gives rise to an invariant subspace 

K{A; x) = {p{A)x\pe V) 

of A, where V denotes the set of polynomials. In the self-adjoint case, the function 
space C(r2) arises again once A is represented on K{A]x) with an antilinear Jacobi 
operator 

where is a complex symmetric (typically infinite) matrix and r denotes the con- 
jugation operator. To characterize the Jacobi parameters on the diagonals of J#, a 
biradially supported theory for C{r2) is devised. A curve in C is said to be bi- 
radial if it intersects every origin centred circle at most at two points. Namely, it is 
shown that to any given set of bounded Jacobi parameters corresponds a compactly 
supported biradial measure yielding these parameters in an iterative orthogonaliza- 
tion process for polynomials p in the respective space. The case of unbounded 
Jacobi parameters can be treated in terms of conditions on the moments recorded in 
a Hankel-like matrix. Regarding the lack of uniqueness of this correspondence, the 
finite dimensional case is completely solved. 

The paper is organized as follows. In Section [2] spectral mapping theorems for 
antilinear operators are derived. The corresponding function space structure is identi- 
fied. In Section [3] a theory for self-adjoint antilinear Jacobi operator is developed. To 
deal with the Jacobi parameter problem, biradial spaces are introduced. In Section 
mthe unbounded case is considered. In Section [5] the non-injective determination of 
the Jacobi parameters is solved in finite dimensional cases. 

2. Functions of antilinear operators and spectral theory. A continuous 
additiv^ll operator on a complex Hilbert space H is real linear, i.e., it commutes 
with real scalars. This fact can be found already in Banach's classic book on linear 
operators [2]I1 Denote the family of such operators by B{H). The norm of i? G B{H) 
is defined as 

||B|| = sup ||i?a;||. 

|x|| = l 

The adjoint of B is the real linear operator B* satisfying 

{Bx,y)^ = {x,B*y)^ 

for every x,y ^ H. Here (•,-)k = Re (•,•), where (•,•) denotes the inner product on 
H . If B* = i?, then B is said to be self-adjoint. 

There exists a unique separation of B into its complex linear and antilinear parts 

as 

B = C + A, 

^An operator S on is additive if B{x + y) = Bx + By for any x,y £ H. 

^It is instructive to bear in mind that Banach's book started with additive operators and dealt 
only with real scalars, a fact which sometimes was regarded as curious |19l p. 397]. 
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where 

C =^{B- iBi) and A = ^{B + iBi). (2.1) 
That is, then C is complex hnear while A is antihnear, i.e., AX — XA holds for any 

A e C. 

Of course, the case of complex linear operators, i.e., when A = 0, has been 
extensively studied. In what follows, we are more interested in the antihnear case, 
i.e., when C = 0. To this end, the following space of polynomials was introduced in 

Definition 2.1. Polynomials of the form 

Hi 

p(A)=^(a2fc + a2fc+iA)|A|2'= (2.2) 

fc=0 

with ttfe G C and, for j even aj+i ~ 0, are denoted byVj{r2). Their union U°°^QVj{r2) 
is denoted by V{r2). 

For a nonzero polynomial, the greatest integer I such that a; 7^ in (|2.2p is called 
the degree of the polynomial. 

To see how these relate with real linear operators, take a standard analytic poly- 
nomial p{X) — X]fc=o "fcA*^- If B e S(i?), then p{B) is defined as 

p{B)^Y.^uB\ 

Of course, polynomials in complex linear operators appear regularly. Polynomials 
in antihnear operators have recently proved useful in numerical linear algebra and 
approximation theory [51[TT1[T^. A first notable difference between these two extremes 
is the fact that a polynomial in an antihnear operator typically becomes genuinely a 
real linear operator, i.e., it has nonzero complex and antihnear parts. Second, Vj{r2) 
is a natural function space to analyze polynomials in antihnear operators. 

Namely, associate with p the polynomial p defined in (|2.2p . At the most fun- 
damental level, such a transformation occurs with a spectral mapping theorem for 
antihnear operators. This is readily seen in terms of eigenvalues. If A is antihnear 
and Ax = Xx for a nonzero x G H and A G C, then (|1.1|) holds. Consequently, the 
eigenvalues of A are mapped with p to be among the eigenvalues of p{A). 

Definition 12.11 was not aimed at maximal generality, however. First, aside from 
eigenvalues, we have a more general notion of spectrum. Second, Laurent series in a 
complex linear operator can be defined as long as the spectrum is contained in the 
annulus of convergence. The spectrum of a real linear operator is defined in a natural 
way as follows. 

Definition 2.2. The spectrum of B G B{H) consists of those points A G C for 
which XI — B is not boundedly invertible. The set of these points is denoted by o^{B). 

The spectrum of a real linear operator is always compact. For antihnear operators, 
it is circularly symmetric with respect to the origin, and it can be empty as well; see 
[13] . Hence, a lack of spectral radius means that spectrum alone is not sufficient to 
determine convergence of Laurent series in real linear operators. Still, if / is analytic 
in an annulus centred at the origin with the Laurent series 
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then, under appropriate assumptions on an antihnear A, the real linear operator f{A) 
is well defined in terms of the seriesH (If — for k — —1, —2, . . ., then we have 
a disc centred at the origin.) We assume that the spectrum (j{A) is contained in the 
annulus of convergence. By inspecting eigenvalues, / gets transformed in the process 
completely analogously as 

CO 

/(A)= J2 i'^^'^+^^k+iX)\M^'=ui\\\') + v{\\\')X. (2.4) 

k— — oo 

Clearly, u and v are polynomials if and only if / is. Observe that / maps origin 
centred circles to circles. In particular, we regard these functions as biradial for the 
following reason. 

Example 1 . The functions u and v in (|2.4p are "biradially" uniquely determined 
in the following sense. Suppose ^1,6*2 G [0, 27r) and 9i ^62- Then, for r fixed, the 
conditions 

J u{r'^) + v{r'^)re^^ = a{r) , , 

\ u(r2) + t,(r2)re«2 = b{r) ^ ' 

with a{r) and b{r) given, determine the values of u(r^) and w(r^)lf| Of course, if w = 0, 
then we are dealing a standard radial function. 

Denote by C(r2) functions of the form on right-hand side of (|2.4p . Assuming that 
the antilinear A is additionally self-adjoint, we shall allow continuous functions of this 
type by noting that A^ is complex linear positive semidefinite. To ensure that this is 
consistent, let us invoke the following proposition. 

Proposition 2.3 CTS"). Let A e B{H) he antilinear. Then A G cf{A) if and only 

z/|A|2efT(A2). 

Definition 2.4. Let A e B{H) he antilinear and self-adjoint. Let u and v be 
complex valued continuous functions defined on the compact subset cr(yl^) C [0, cx)). 
Then the function f : cr{A) — > C defined by 

/(A)=u(|An+^;(|A^A (2.6) 

is called a continuous hiradial function. Define 

f{A)=uiA') + v{A')A. 



Clearly, C(r2) is a vector space over C. There exists a natural notion of (non- 
commutative) product as well. To this end, consider again (II. ip . For two elements 
/(A) = wi(|A|2) -f wi(|A|2)A and 5(A) = uadAp) + t;2(|A|2)A of C(r2), we have for the 
eigenvalues 

f{A)g{A)x = h{X)x = (/ * g){X)x (2.7) 

once we define 

MA) = um>2i\\\')+\\fvii\x\')Mm+ {um')M\M')+^Mmv2{\\\')) a. 

•^Form limsupj._^_|_^ and compare with Iimsupj._^_|_^ ||-'^/-'. 

*This makes extending C(r2) apparent. To determine tlie values 2fc-radially, we are lead to 
consider functions of the form X]j=o^ '^j (l-^P*)-^"' 1 where Uj are sufficiently smooth and k £ N. 
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Or, in short 

f*g:^{ui+viXT){u2 + V2X), (2.8) 

where r denotes the conjugation operator. 

After identifying an appropriate function space for antilinear operators, we are 
concerned with having a spectral mapping theorem. To this end we need the fohowing 
lemma. 

Lemma 2.5. Let A e B{H) be antilinear and self-adjoint, and let f : <7{A) — > C 
be a continuous biradial function. Then for any e > there exists a polynomial p such 
that 

\\f{A)^p{A)\\ < e and max |/(A) -p(A)| < e. 



Proof. The result follows by applying the Weierstrass approximation theorem to 
u and V in (12.61) on the compact set cr(A^) and using the continuous function theory 
for the complex linear A^. □ 

Theorem 2.6. Let A G B{H) be antilinear. Assuming f{A) is defined for f in 
(|2.3p . there holds 

f{a{A)) C a{f{A)). 

Moreover, if A is self-adjoint, then this further holds for all continuous biradial func- 
tions f : a{A) -^C in JM]). 

Proof For any B e B{H) holds 

aaiB)U^jB^^a{B), (2.9) 

where aa{B) denotes the approximate point spectrum of B and the bar denotes com- 
plex conjugation (not closure). 

Take A E a-a{A). Then for any given e > there exists a unit vector x E H such 
that Ax ^ Xx + v with ||u|| < e. Let € R. Since Ae^^x = e"*^Aa; + e~^^v, we may 
conclude that aa{A) is circularly symmetric with respect to the origin. Therefore by 
()2.9p . also (J (A) is circularly symmetric with respect to the origin. 

Consider a rational function r(A) = J2k=-i '^kX^ with l,j G N. (For this the iden- 
tity (|l.ip holds analogously.) By continuity, we have r{A)x — f{X)x -\-w, where 
can be made arbitrarily small by choosing e small enough. Consequently, f(A) G 
aa{r{A)) whenever A € cra{A). Regarding /, we have ||/(^) — ''(^)|| < e and 
max^go-(yi) |/(m) ~ ^(a*)I ^ £ foi' any e > by choosing I and j large enough. Thereby 

\\f{A)x - f{X)x\\ < \\f{A)x - r{A)x\\ + \\r{A)x - f(A)a;|| + ||f(A)x - /(A)x|| 
< 2e+ 

and it follows that /(A) S (Ja{f{A)). The case with self-adjoint A is similar, only then 
we have a polynomial in place of r obtained by applying Lemma 12.51 

Similarly can be dealt with the case A S aa{A*). To see this, it suffices to consider 
the case of an eigenvalue of A* and the corresponding circle of radius |A| centred 

at the origin. The claim follows from r{A)* — J2i=-i^*'^'^ = 1 1 (^^2! + 

a2k+iA*)A*'^'^ after multiplying by the corresponding eigenvector. □ 
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In general, the equality does not hold. For example, if the spectrum of A is 
empty, then with /(A) — |Ap we have f{a{A)) = 0, but cr(A^) ^ 0. The equality can 
be established in the following case, whose proof will be postponed until the end of 
Sectioning 

Theorem 2.7. Assume A G B{H) is antilinear and self-adjoint. Then, for all 
continuous biradial functions f : cr{A) — !■ C m (|2.6|) . 

f{a{A))=a{f{A)). (2.10) 

Corollary 2.8. Assume A G B{H) is antilinear and self-adjoint. Then 

||/(A)|| = max |/(A)| = max |A|. 
\e>j{A) ag<t(/(A)) 



For an illustration, Hankel operators yield an immediate nontrivial family of self- 
adjoint antilinear operators. 

Example 2. Hankel operators [17] constitute a natural family of self- adjoint 
antilinear operators once treated as follows. (See also [HI Sec. 7 and 8].) Denote by 
-ff^(D) the Hardy space, where D is the unit disc. Let a € L°°{T), where T denotes 
the unit circle and P : L^(T) _ff^(D) is the orthogonal projector onto _ff^(D)l£| 
Then define 

g ^ PMag 

on 77^(D), where Ma is the multiplication operator g i— > ag. Represented on we 
have Ht, where H is an infinite Hankel matrixlf] 

Like Hankel operators, self-adjoint antilinear operators have also been studied in 
a complex linear setting more generally; see [7j[8] where also many examples are given. 

We do not know how to completely characterize those antilinear operators for 
which the spectral mapping Theorem l2. 61 holds with equality. However, we conjecture 
that these are precisely those antilinear A for which cr(yl^) C [0, cx)). Necessity follows 
by taking /(A) = A^ and using Proposition 12.31 In finite dimensional H, sufficiency 
is established by the notion of contriangularizability, which is equivalent to tT(A^) C 
[0,oo) 9 . (For the probability of being contriangularizable, see [12].) 

Example 3. Suppose an antilinear A € B{H) is such that ct{A^) C [0,oo). 
Then also Gelfand's formula 

lim ll^-'ll^/^ = max |A| 

holds. 

3. Antilinear Jacobi operators and theory for C(r2). Like in the com- 
plex linear case, the spectrum is intimately related with the notion of invariant sub- 
space. That is, if an antilinear operator has an eigenvector, then its span yields an 
invariant subspace. By invariance is meant the following. 

Definition 3.1. A subspace K of H is said to be invariant for an operator 
B e B{H) if BK C K. 



^Recall that ^^(D) can be identified with those elements of (T) which have vanishing Fourier 
coefficients for negative indices. 

^This antilinear treatment leads to different problems. For instance, unitary similarity for Ht 
becomes unitary consimilarity for H . We are only aware of the unitary similarity problem for H I14| . 
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For an antilinear A e 13(H) and x £ H, form the invariant subspace 

KiA;x)^{piA)x\peV}, (3.1) 

where V denotes the set of polynomials. If x is an eigenvector, then K(A; x) is one 
dimensional. For the other extreme, if K{A; x) — H, then x is said to be a cyclic 
vector for A, like in the complex linear case. 

An orthonormal basis oi K{A; x) can be generated analogously to the finite dimen- 
sional case described in [5]. If the dimension of is infinite, then A is represented 
as 

H^T : f{N) /2(N), 

where the (infinite) matrix is of Hessenberg type with real subdiagonal entries 
and T denotes the standard conjugation operation on /^(N). If the dimension of p.ip 
is finite, say n, replace P{N) with C". In case A is self-adjoint, is a tridiago- 
nal complex symmetric Jacobi matrix. In this way complex Jacobi operators arise 
naturally in antilinear operator theory. 

Proposition 3.2. Let A e B{H) be antilinear and self-adjoint, and let x,y & H. 
IfyJ.K{A;x), then K {A; y) ± K (A; x) . 

Proof. Let p and q be polynomials, and denote q{X) ~ u(A^) + w(A^)A, where u 
and V are polynomials. Then 

ipiA)x,q{A)y) ^ {p{A)x,u{A^)y) + {p{A)x,viA^)Ay) 

= iu{A^)p{A)x,y) + {Av{A^)piA)x,y) 

= (p(A)x, y) + {q{A)x,y) = 0, 

where we have used the fact that is C-linear self-adjoint and p, q are some poly- 
nomials. □ 

By Proposition 13.21 given a bounded self-adjoint antilinear operator A on H, we 
can express H as an orthogonal direct sum of invariant subspaces 

i/ = 0i^(A;a;„). (3.2) 

a 

We denote the restriction of A to these subspaces by 

— (3.3) 

We have 

a{A) = \Ja{A^), (3.4) 



since <j{A'^) — [J^ a{A'^), and we apply Proposition 12.31 noting that <t{A) and a{Aa) 
are circularly symmetric with respect to the origin. 

With these preliminaries, next we show that there exists a correspondence between 
self-adjoint antilinear Jacobi operators and positive measures on the plane, much as 
in complex linear operator theory. The construction relies on theory for C(r2) 
by orthogonalizing polynomials (12.21) supported on these measure spaces. Since the 
invariant subspaces K{A; x) may be finite or infinite dimensional, we shall deal with 
them separately, starting with the finite dimensional case. Then, once we have handled 
the infinite dimensional case. Theorem 12 . 71 will be proved. 

It is instructive to bear in mind that, as described in Example [3J everything 
applies to Hankel operators. 



8 



M. HUHTANEN AND A. PERAMAKI 



3.1. Biradial measures with finite support and antilinear Jacobi opera- 
tors on finite dimensional spaces. In what follows it is shown, by partly following 
[T5), that there exists a correspondence between (discrete) biradial measures and an- 
tilinear Jacobi operators. For convenience, the following notation for the monomials 
appearing in Definition 12.11 is employed. 

Definition 3.3. Let k €N. Denote hy{\)^ the monomials defined by 




A|Ap^ , if k is odd and k = 2j + 1, 
|A|'', if k is even. 



To put it short, biradiality means the second item in the following definition. 
Definition 3.4. Assume n G N. Lei pfc > 0, Aj. e C (fc = 1, . . . ,n) be such 

that 

(i) YTk=iPk = 1, 

(ii) Ai, . . . , A„ are distinct and any origin centred circle intersects at most two of 
them. 

Furthermore, the complex numbers Ai, . . . , A„ are assumed to be ordered such that 

|A2fe-i| = |A2A;| fork = l,...,m, 
|A2fe-i| < |A2fc+i| fork = l,...,m- 1, 

|A2m+l| < |A2m+2| < • • • < |A,i| 

for m e N. Then a positive measure on the plane defined by 

n 

P = '^PkSx^ (3.5) 

fc=i 

is called a biradial measure with finite support. 

Assume p is a biradial measure with finite support. On V{r2) let us use the 
inner product 



/ pqdp = y p{\k)q{Xk)pk. (3.6) 



k=l 



Consider the monomials 1, (A)^, (A)^, . . . , (A)" \ Executing the Gram-Schmidt or- 
thogonalization process with respect to this the inner product yields an orthonormal 
sequence of polynomials po,Pi, ■ . . ,Pn~i such that each pj has degree j. This is a 
consequence of the following proposition guaranteeing that = > for all 

p S V{r2) with degree less than n. 

Proposition 3.5 (T^). Letp e Vd{r2) be nonzero. The following claims hold: 

1. If p has two distinct zeroes of the same modulus, then all numbers of that 
modulus are zeroes. 

2. Let m be the number of nonzero moduli for which all numbers of that modulus 
are zeroes and let s be the number of moduli for which exactly one number is 
a zero. Then 2m + s < d. 

Consider the product (|2.8p . Expressing the polynomial Arpj(A) — Apj(A) as a 
linear combination of pq, . . . ,Pj+i by imposing orthogonality gives rise to the three 
term recurrence 



/3j+iPj+i(A) = Apj(A) - aj+iPj{X) - /3jPj_i(A), (j = 0, . . . , n - 2), (3.7) 
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where p_i(A) = 0, ctj+i = {^PJ^Pj) — i^pJ^Pj+i) ■ Observe that > 

since the leading term of Pj+i has positive coefficient, a consequence of executing the 
Gram-Schmidt process. Proposition 1331 the space {Vn{r2), (•, •)) is n dimensional. 
Therefore Ap„_i(A) is a linear combination of po, . . . ,Pn~i and hence /?„ = 0. In this 
way we have associated these so-called called Jacobi parameters ai, . . . , a„ € C and 
/3i, . . . , Pn-i > with the given biradial measure p. 

To express this linear algebraically, let Q — (qkj) G C"^" be the unitary matrix 
with columns j defined by 



qkj = y/pkPj-iiXk), 
Then, by using the recursion p.7p . we get 



{k= l,...,n). 



(3.8) 
(3.9) 



where 



Oi\ Pi 
Pi "2 

... 
... 

£)# = diag(Ai, . . 





, A„). 






ttn-l Pn~l 
Pn-1 OLn 



(3.10) 



(3.11) 



In this way, given a biradial measure with finite support, we have the corresponding 
complex Jacobi matrix (|3.10|) . 

Suppose, conversely, that we are given a complex Jacobi matrix p.lOp and we 
want to find a corresponding biradial measure. We start with the following proposi- 
tion. 

Proposition 3.6. Let ai, . . . , a„ e C and let . . . , Pn-i > 0. Then there exist 
numbers Ai, . . . , A„ G C and a unitary matrix Q G C"^" with positive first column 
such that the equation (|3.9|) is satisfied, when and are defined by formulae 
p.lOp and p. lip . Furthermore, the numbers Ai,...,A„ satisfy the property (ii) in 
Definition \S.4\ 

Proof. By f9l Corollary 4.4.4], there exists a unitary matrix U E cix" and a 
nonnegative diagonal matrix E# — diag(tTi, . . . , (t„) such that 



(3.12) 



We show that no three of (Ti , . . . , (T„ can be equal. To reach a contradiction, we 
assume, without loss of generality, that ci = cr2 = fa- Let V G C^x" be the first 
three rows of U and let vi,. . . ,Vn be the columns of V. By (|3.12p . we have the 
recursion 



aivi = [3iV2 + aivi, 



1, 



, ,n- 2. 



It follows that each Vj is a linear combination of vi and vi . This implies that the rank 
of V is at most two, contradicting the fact that the rank is three. 



Let e G 



be a diagonal matrix such that Q 



U has a nonnegative first 



column. Define = e^'®S#. Then the identity (|3^ holds. That ah Ai,...,A„ 
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are distinct in p.lip can be shown by an argument similar to the one given above. 
Furthermore, if qj is a column of Q*, then J^TJJ ~ XjQj and it is easy to see that the 
first entry of qj must be nonzero. □ 

Corollary 3.7. Let ai, . . . ,a„ e C and let Pi, . . . ,/3„_i > 0. Then there exists 
a biradial measure fi with support of cardinality n such that the complex Jacobi matrix 
corresponding to fi is given by (|3.10|) . 

Proof. Let Q = (qkj) € C"^" and e C"^" be according to the previous 
proposition and let pk = q\i. Define the measure p = Y^^^i Pk^Xk- For each j = 
1, . . . ,n we then regard p.8p as a system of equations with Pj-i as the unknown. 
These systems have invertible Vandermonde-type matrices, by Proposition 13. 5[ and 
are therefore uniquely solvable. A glance at the equation (13.91) then shows that the 
polynomials Pj-i satisfy the recursion p.7|) with the same coefficients as the given 
complex Jacobi matrix J^. □ 

To sum up, based on constructing an theory for C(r2) and then orthogonalizing 
polynomials (j2.2l) according to p.7p . we have shown that the mapping from biradial 
measures with finite support to complex Jacobi matrices of type given in Corollarv l3.7l 
is surjective. It is not injective, however. This lack of injectivity will be completely 
described in Section [Sj after establishing surjectivity in the infinite dimensional case. 

3.2. Biradial measures with infinite support and antilinear Jacobi op- 
erators on infinite dimensional spaces. We now turn to biradial measures with 
infinite supports in preparation to the study of antilinear Jacobi operators on infinite 
dimensional Hilbert spaces. We denote R"*" = [0,oo). 

Definition 3.8. A Borel probability measure p on C is said to have finite mo- 
ments, if it satisfies 

[ \X\" dp{X) < oo for all n = Q, 1,2,.... (3.13) 
Jc 

Definition 3.9. Let Pi,P2 be positive Borel measures on R+ such that 

p+(R+)+p+(M+) = l. 

Let <pi, <p2 be real-valued Borel-measurable functions on and define the measurable 
transformation R^. : C — C by 

Define the positive Borel measures pj on C by 

p,{E) ^ p+{R^^ {E) n (K+ X {0})), (j = 1, 2) 

where E d C is Borel-measurable. Then p = pi + P2 is called a biradial measure. 
Moreover, it is called symmetric if <j)i{r) — (j)2{'r) + tt (modulo 2'k). 

We can represent a biradial measure p in an equivalent form as follows. Let p+ = 
pt + Pt ^^'^ ^1,02 be the Radon-Nikodym derivatives aj = dp'j' /dp~^ (j = 1,2). 
Note that < ai{r),a2{r) < 1 and ai{r) + a2{r) = 1 for almost every r. Suppose 
i? C C is measurable and write as follows (where we use R^. (E) as short-hand for 
R^.{E) n (M+ X {0}) with the vacuous second coordinate removed) 

/>00 pQO 

p,{E)^ XR,,(E)dp+= a,{r)6r{R^^{E))dp+{r) 
Jo Jo 




a,{r)5^^.,,^.,{E)dp+{r). (j = l,2) 
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Hence we can write 

POO 

p{E)= ^ir{E)dp+{r), where = V a,(r)5^^,,,M (i?). (3.14) 

The formula p.l4p is the disintegration of the measure p with respect to the measur- 
able map T : C — ?► M+, T{z) = \z\, and [4'. Conversely, we could define a measure 
by the formula p. 141) and rewrite it in the form of Definition 13.91 

Example 4. Let /i be a positive Borel measure on K. Define the measures 
p+ on R+ by p+{E) = p{E) and p+{E) = n{-E D (-oo, 0)), where E C M+. Let 
= and 02 = T^- Then the measure p on C in Definition 13.91 is supported inside 
R X {0}, and p{F x {0}) = n{F) for all Borel sets F C M. 

Example 5 . Let pf be a Borel probability measure with support R+ and abso- 
lutely continuous with respect to the Lebesgue measure. Let {^fc}^i be a partition 
of R"*" into metrically dense subsets, i.e., for all k we have Ak Borel measurable and 
ptiA-k n {x - e, X + e)) > for all x e R+, e > (e.g. [6]). Let {(?fc}^i be an enumer- 
ation of the rationals and define 0i = X^fc^i Ik^A^- Let — 0. Then the measure p 
in Definition 13.91 satisfies p{B{z,r)) > for every open disc B{z,r) C C and therefore 
the support of p is C. 

With respect to the monomials (A)*^, measures on C can be reduced in dimension 
as follows. 

Lemma 3.10. Let p, he a Borel (probability) measure on C with finite moments. 
Then there exists a symmetric biradial measure p such that 

pdp= pdp for allpe V{r2). (3.15) 
Jc Jc 

Moreover, if p is compactly supported, then also p is. 

Proof. Define the measurable map T : C R+, T{z) = \z\ and let {pr}r£R+ be 
the disintegration of p with respect to T and the image measure = pT^^ 4 . For 
nonnegative integer fc, we then have 

/ (A)%(A)= r f {Xfdpr{X)dp+{r), 

Jc Jo J\\\=r 

where, for almost all r, pr is a Borel probability measure supported in {|A| — r}. For 
odd k, write k ^ 21 + 1, and then we get 

/(A)"dp(A)=/ / Xdpr{X)r^'dp+{r), 

Jc Jo J\\\=r 

and we denote the inner integral (center of mass) on the right-hand side by C(r). 
Then we can choose oi (r ) , 02 (r) such that < ai(r),a2(r) < 1 and ai(r)+a2(r) = l, 
together with (j)i{r),(f>2{r), such that 

2 

C(r) = ^aj(r)re*'^^W. 
i=i 

If C(r) 7^ is on the circle {|A| ~ r}, we can set 02 (r) — and then 0i(r) is uniquely 
determined (modulo 27r). If C{r) 7^ is inside the circle, these choices are non-unique 
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as illustrated in Figure [??Tl and we let 02 (^) — 0i(^)+^ (modulo 27r) to get symmetry. 
If C(r) = 0, we let (/'i(r) — 0, 02 (^) = tJ"- We now define the measures 

2 

(Note that if C(r) 7^ 0, our choices for ai(r), 02 (r), 0i(r) and 02 (f") make the unique 
measure of this form.) We observe that 



AdMr(A) = C(r) = / \dVr{X). 

|A|=r J|A|=r 

The measure defined by p{Ej) = i'r{E) dp^{r) is symmetric biradial, see formula 
p.l4p . and we now have 

\\f dii{\) = j {Xfdp{X) (3.16) 



for all odd k. It is easy to see that p. 161) is true for even k as well, and then (|3.15p 
follows. 

The final claim follows from the fact that is compactly supported if is. □ 

Proposition 3.11. Let p, be a compactly supported biradial measure. Then the 
set of polynomials V{r2) is dense in L^(/i). 

Proof. Let i? > be such that the support of p is contained in an origin-centred 
closed disc of radius R. We denote the disintegration of p as in formula (|3.14p , and 
denote Ai(r) = re"^^^''\ \2{r) = re'-^^^''). 

Take / S L'^{p) and e > 0. Let g be a compactly supported smooth function on 
C such that ||/ - 5||^2(^) < e/2. For r such that Ai(r) = A2(r), let ^(r^) = g{\i{r)) 
and v{r'^) — 0, and otherwise let u{r'^) and w(r^) be the unique numbers such that 

g{\,{r))^u{r^)+v{r^)\,{r), (j = l,2). 

For r such that Ai(r) 7^ A2(r), we now have 

ff(Ai(r))-5(A2(r)) 



Ai(r)-A2(r) 



= v{r'), 



where the left hand side is bounded. Hence u is a bounded function and it follows 
that u is bounded as well. Let p and q be ordinary polynomials such that 



\u{r)-p{r)\dp'^{r)<- and / \v{r)r-q{r)r\dp'^{r)<-. 

8 " 



Then we have 

' \f{\)-p{\\?)-q{\\?)\? dp{X)<e. 



□ 

With these measure theoretic preparations, let p be a Borel (probability) measure 
on C with finite moments. On V{r2), we define the inner product 

{p,q) = I pqdp. (3.17) 

Jc 
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We then apply the Gram-Schmidt process to the monomials 1, (A)^, (A)^, . . . with 
respect to the inner product p.6p and obtain an orthonormal sequence of polyno- 
mials po,pi,p2, ■ ■ ■ , where pj has degree j. The process breaks down if and only if 
dim7'„(r2) < n for some n. In the breakdown case, for the least such n, we have 
the orthonormal polynomials po, ■ ■ ■ ,Pn-i and the corresponding Jacobi parameters 
{Q!j}"^i, similar to the case of biradial measures with finite support. If the 

process does not break down, we get infinitely many Jacobi parameters {aj}j°^i with 

positive {/3j}j*Li and the recursion p. 71) holds for all j = 0,1,2, These Jacobi 

parameters are recorded in the infinite matrix 



ai f3i 
Pi a2 h 
/32 as 



(3.18) 



Proposition 3.12. Let p,p' he Borel (probability) measures with finite moments 
and at least n Jacobi parameters, {oLj}"^^^, {PjYj^l and {a^}"^!, ; respec- 

tively. Then aj — a'j (j = 1, . . . ,n) and l3j = P'j {j — 1, . . . ,n — 1) if and only 

{Xf dp = {X}'" dp' for all k = 0,1,..., 2n-l. (3.19) 
Jc Jc 



Proof. The proof for the real Jacobi parameter case carries over easily [^TJ Propo- 
sition 1.3.4]. □ 

The following is a version of Favard's theorem for bounded antilinear Jacobi op- 
erators. We prove it similarly as in |21) . 

Theorem 3.13. Let {aj}'^i, {j3j}J^i be bounded Jacobi parameters. Then there 
exists a compactly supported symmetric biradial measure p on C .such that the Jacobi 
parameters corresponding to p are {aj}'j^i, {/3j}^i- 

Proof. By Corollary 13.71 for each n there exists a biradial measure /i„ such 
that the corresponding Jacobi parameters are , and we denote the 

In] . 

corresponding complex Jacobi matrix by J# . Since the n points supporting pn are 
coneigenvalues of J^\ the support of pn is contained in the closed disc B{0, || j]r'||) C 
5(0, II J#r|| ) =: B. Moreover, Pn{C) = 1 is bounded for all n and therefore there 
exists a weakly converging subsequence {pn^yfLi such that 

lim / fdpn^ ^ i fdp 
j-^ooJb Jb 

for all continuous functions / on B, where p is a Borel probability measure supported 
in B. For each nonncgative integer k, we choose /(A) = (A)*^ and use Proposition 
l3.12l to conclude that the Jacobi parameters of p are {ajjj^i, Moreover, an 

application of Lemma 13.101 replaces p with a symmetric biradial measure. □ 

We next obtain a spectral theorem for bounded antilinear self-adjoint operators. 
For an approach using spectral integrals, see [20]. 

Theorem 3.14. Let A be a bounded antilinear self-adjoint operator on a Hilbert 
space H . Suppose there exists a cyclic vector of A. Then there exists a compactly 
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supported symmetric biradial measure n on C, and a C-linear isometric isomorphism 
U : H ^ L'^{^), such that 

UAU-^ = At, (3.20) 

where At denotes the multiplication operator /(A) A/(A) on L^(/i). 

Proof. Let x £ H he the cyclic vector and let {qj}J^i be the orthonormal basis 
of H generated by the Arnoldi process starting from the vector x. Moreover, let 
be the generated complex (infinite) Jacobi matrix with parameters {aj}J^i, {I3j}°°^i. 
By Theorem 13. 131 the measure fi exists, with orthonormal polynomials {pj}'j°^Q, such 
that 

PjPj-iW + (^j+iPjW + /3j+iPi+i(A) = Apj(A) (j 0, 1, 2, . . . ), 

where p_i = 0. The set of polynomials {pjJj^Q is an orthonormal basis of L^ifJ-) 
by Proposition 13. 1 ll We define the isometric isomorphism U by U{qj) ~ Pj-i (j = 
1,2,...) and the claim follows. □ 

Corollary 3.15. Under the assumptions of Theorem \3.14\ for any continuous 
biradial function f defined on cr(^), 

Uf{A)U-'^u{\X\') + v{\X\')XT, 

where u and v are as in (|2.6|) . 
Proof. Apply Lemma 12.51 □ 

Proposition 3.16. Biradial measures fi given by Theorem \3. 14\ are supported in 
the spectrum cr(A). 

Proof. By Theorem 13. 14[ we have UA'^U^^ = l-^P and therefore 

Up{A^)U-' ^ p{\X\') (3.21) 

for all ordinary complex analytic polynomials p. The operator A^ is C-linear self- 
adjoint and positive semidefinite. 

Let Ao ^ cr{A). Let g{X) — /(|Ap), where / is a nonnegative compactly supported 
function on such that / vanishes on cr(A^) and (?(Ao) — 1, which is possible due to 
Proposition [231 Let {pj} be a sequence of ordinary polynomials converging uniformly 
to / on the support of /. Then \\f{A'^)—pj{A^)\\ ^ as j -^• cx), and, since f{A^) = 0, 
from p.2ip we have 

/ g(A)dAi(A)=0. 

Hence there exists an open set V C C containing Ao such that fJ.{V) = and therefore 
Ao is not in the support of /i. □ 

We are now ready to prove Theorem 12.71 and Corollary 12.81 

Proof of Theorem \K7\ The inclusion C for (|2.10p is Theorem [221 For the 
inclusion D, take Aq ^ /((t(A)), and let H and A be decomposed as in <\3.2\i and (|3.3p . 
By Corollarv l3.15[ we have 

;7(Ao/a - f{Aa.))U-' = Ao - w(|An - v{\X\'')Xt. (3.22) 

Since cr(A) is a compact set. 



c :— min 

xe<j{A) 



= min |Ao-/(A)|>0, 
Aec7(A) 
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and, since f{(j{A)) D f{a{Aa)) by ([5^ . we have 

0<c< min \Xo~u(\X\'^)\~\v(\X\^)X\ 

AG<t(Aq) ' III 



for all a. 



Hence the multiplication operator on the right-hand side of equation p.22p is bound- 
edly invertible, with the inverse having the norm estimate 



max 



\xo~u{m\-\v{\xm 



1 

< - < oo 
c 



for all a. 



We have shown that the operators XqIc — /(^q) are boundedly invertible for all a, 
with the norms of their inverses uniformly bounded. Hence XqI — f{A) is boundedly 
invertible, so Xq ^ a{f{A)). □ 

Proof, of Corollary \2.8[ The second equality follows from Theorem 12.71 For 
the first, let H and A be decomposed as in p.2p and p.3|) . and note that ||/(^)|| = 
sup„ \\f{Aa)\\. Then, by Corollary [HSl 

< / {\u{\X\')\ + \vi\Xmr\h{X)\' df,{X) 

< sup |/(A)p||;i||2. 

where we denoted h = Ug. Hence < sup;^^^^^^-) l/(-^)l- The inequality 

in the other direction is established by considering invertibility of XIa ~ f{Aa) = 
A(/q — X~^ f (Aa)) via the Neumann series expansion in the usual way. □ 

Starting from [3], multiplication operators have played a major role in under- 
standing the C-linear invariant subspace problem. For the corresponding problems 
in the antilinear case, the role of multiplication is taken by the operator p. 201) . 

4. Unbounded antilinear Jacobi operators and the moment problem. 

Let p be a biradial measure with finite moments and define 



TOfe = / (A) dp{X), 



fc = 0,1,2,. 



(4.1) 



In the corresponding moment problem we are given an arbitrary sequence of complex 
numbers {wfe}^^ such that, after a possible scaling, mo ~ 1. Of course, then neces- 
sarily m2fc > for k € N. The problem consists of finding a biradial measure p such 
that the equation (|4.ip is satisfied. The answer can be given in terms of the matrix 



M 



Too TTll 
TOl TO2 
TO2 TOj 



TO2 TO3 
TO3 TO4 
TO4 



TO3 TO4 TO5 TOg 



(4.2) 



Observe 



where, denoting again the conjugation operator by t. Mi. 
that AI is Hermitian, as well as, in a certain sense, a Hankel-like matrix. Thereby it 
is natural to contrast this with the classical Hamburger moment problem. 

Theorem 4.1. Let {mk}'^i C C and toq — 1. Then there exists a biradial 
measure p with finite moments such that (|4.1I) holds if and only if the suhmatrix 
Mi;k,i:k of (|4.2p is positive semi-definite for all k. 
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Proof. That the claim holds with a Borel probability measure fj,, in lieu of a 
biradial measure p, follows from 22, Theorem 1]. This measure p can then be replaced 
by a biradial measure p by Lemma 13.101 □ 

Corollary 4.2. The mapping from biradial measures to complex Jacobi param- 
eters is surjective. 

Proof. Let {ctj, f3j}J^i G (C x (0,oo))°° be complex Jacobi parameters. Define a 
sequence of polynomials {pj}'j°^o by the recursion (|3.7p starting with pq = 1. Define 
an inner product in 'P{r2) by 

{p,q)=J2<^jb], (p,geP(r2)), 

j>0 

where p — J2j '^jPj ^-i^d q — J2j ^jPj i ^^-d let 

mk = ({X)\l) for fc = 0,1,2,.... 

From the recursion p.7|) we see that multiplication by Ar in V{r2) is represented by 
J^T in the basis {pj j^oi where is the infinite complex Jacobi matrix given by 
formula (|3.18p . From this we find 

m,+, = eJ{J^ry+=e, = r= ((J*r)^"ei)* ((J*T)'ei) = r= ((A)\(A)^) , 

where ei = [l ■ ■ 'j^- Since the inner product is positive-definite, the condi- 
tion of Theorem 14.11 is satisfied and therefore there exists a biradial measure p such 
that 



{^(A)^ 1^ = mfe = J {\)'"dp for all k. 



Since r^ (A)'^"' = (A)'(A)"', we get {p,q) — J^pqdp for any polynomials p and q. This 
completes the proof. □ 

5. Noninjectivity of the mapping of measures to Jacobi operators. The 

mapping of biradial measures to antilinear Jacobi operators is surjective, but non- 
injective (even for measures with bounded support). We do not discuss this in full 
generality, but content ourselves with a precise characterization in the finite dimen- 
sional case covered in Section [3?T] 

Theorem 5.1. Let p,p' be biradial measures with supports of cardinality n. 
Then p and p' have the same Jacobi parameters if and only if (using the notation of 
Definition \3.4^ we have 

m ~ m', 

Pk = p'ki when 2m + 1 < k < n, 

P2k-i + P2k = P2k-i + P2fc: "^^en 1 <k <m, (5.1) 



Afe — A'j,, when 2m + 1 < A; < 

l-^2fe| = lAafel, whenl<k<m, 
P2fe-iA2fc-i + P2kMk = P2fe-i-^2fc-i + P2k^2k^ "^^en 1 <k <m. 
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Fig. 5.1. An illustration of the case where the centre of mass C resides inside the circle. Here 
Pi, p'j^ are the masses of the points and A^, are their positions. The total mass and the centre of 
mass are unique. 

Proof. It is easy to verify that the conditions (|5.ip imply p. 191) . Hence sufficiency 
fohows from Proposition 13. 121 

To prove necessity, assume the condition p.lQp holds. First, the equation p.9p 
implies that J# J# and D^D^ have the same eigenvalues. Hence the sets {|Ai|}"^i and 
{|Ai|}iLi are the same, including multiplicities, which implies the first and the fifth 
condition in (I5.ip . Moreover, the conditions given by p. 191) for even k then imply the 
second and the third condition in (I5.ip . Finally, for odd k, (|3.19l) imply the fourth 
and the sixth condition in (|5.ip . □ 

The preceding theorem says that the total mass on each origin centred circle is 
unique and the centre of the masses on each such circle is unique as well. If the 
centre of mass is on the circle, then all the mass is concentrated on a single unique 
point. If the centre of mass is inside the circle, then the mass is distributed among 
two non-miique points as illustrated in Figure [5Tl 

There is also an immediate application of this result in numerical analysis. 

Example 6. Assume that Af# e C"^" is a very large matrix and & e C". The 
R-linear GMRES (Generalized Minimal Residual) method is an iterative method for 
solving the linear system 

M^x = b. 

In |12j its convergence behaviour was analyzed. If complex symmetric, then it 
can be assumed that is actually a diagonal matrix. Suppose, moreover, that 
has distinct diagonal entries such that each origin centred circle intersected either two 
of them or none, and that and b E K" is a vector with all its entries ones. Then in 
|12i Section 5] it was observed that the numerical convergence behaviour was such 
that the residual dropped only at every other iteration step. 

We can now explain this observation as follows. The norm of the residual vector 
Tfc = 6 — M^'xk, where Xk is the approximation at the fcth step, satisfies 

llrfcll = min \\p{M^)b\\ < min max \p{X)\ \\b\\. 

p(0) = l p(0) = l 

By Theorem I5.1[ there exists a unitary matrix Q e C"^" such that p = Q*b G M" 
and that the matrix — Q*M^Q is a diagonal matrix with the property that if 
A G C appears as a diagonal entry then also —A does. Now we have 

llrfcll ^ min \\p{D^)p\\ < min max |p(A)| 

peVk{r2) peVk(r2) \ec(D#) 

p(0) = l p(0) = l 
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If p(A) = w(|Ap) +w(|Ap)A is the minimizing polynomial to the problem on the right- 
most side, then due to symmetry, also M(|Ap) — ?;(|Ap)A is. Hence, by uniqueness, 
V = and p{X) = u{\X\'^). Thereby we may expect that the residual drops only at 
every other iteration step. 
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